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ABSTRACT 

We present three spectra of Cygnus X-l in the low/hard state from 3-1000 keV 
produced with RXTE and OSSE. We demonstrate that a pure thermal Comptoniza- 
tion model, with Compton reflection fits the data well, with coronal temperatures of 
about 90 keV, optical depths of about 1.3, and reflection fractions of about 0.25 x 2ir. 
We find that no additional components are necessary to fit the data, although the pres- 
ence of an additional non-thermal component does result in a marginal improvement 
to the goodness of the fit. Two of the observations are fit better by an almost purely 
non-thermal Comptonization model with Ti n j of ~ 3.2 for the electrons and the other 
parameters essentially the same as for the thermal model, while the third observation 
is also consistent with this model. Observations that could break this degeneracy are 
discussed. We also demonstrate that the low reflection parameter (i.e. R/2tt < 1) is 
due to a lack of strong curvature in the 30-50 keV range of the spectrum and not 
due to ionization effects. The spectrum changes very little over the three weeks of the 
observations, in accordance with a three week correlation timescale found in the All 
Sky Monitor data. We show that for a purely thermal corona, the radius of the corona 
must be between 3 x 10 6 and 1.5 x 10 10 cm, while the radius of a hybrid corona or a 
purely non-thermal corona is constrained to be between 10 8 and 10 9 cm. 
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1 INTRODUCTION 

Cygnus X-l is a binary star system consisting of an 09.7 
lab supergiant (Gies & Bolton 1986) and a compact object 
of at least 6 solar masses (Dolan 1992). The compact object 
accretes matter from a focused wind off the supergiant star 
(Gies & Bolton 1986) and is one of the brightest persistent 
sources of X-rays in the sky. Cyg X-l is the first system to 
show strong evidence of containing a black hole and remains 
one of the strongest candidates for being a black hole. 

Phenomenologically, the spectrum of Cygnus X-l above 
~ 3 keV is well described as an exponentially cutoff power 
law with a photon spectral index of F ~ 0.6 — 0.7 and a cutoff 
energy of around 150 keV. Additionally the spectrum shows 
a broad excess around 30-50 keV due to Compton reflection, 
an iron K edge above 7 keV and an iron Ka flourescence 
emission line (see e.g., Ebisawa 1996). The hard power law 
has generally been described as arising from thermal Comp- 
tonization of soft photons in a hot, photon starved corona 



(Sunyaev & Triimper 1979; Gierlinski et al. 1997; Poutanen, 
Krolik, & Ryde 1997; Dove et al. 1998). 

A previous attempt to model the X/7 spectrum of 
Cygnus X-l relied on data from Ginga and OSSE (Gierlinski 
et al., 1997). They showed that a thermal Comptonization 
model with kT ~ 100 keV, r ~ 1 plasma could fit the data, 
but argued that patterns found in the residual plot suggest 
a two component model with a hot moderate optical depth 
(t ~ 1 — 2, kT ~ 100 keV) plasma and an additional Wien- 
like component with kT ~ 50 keV and r ~ 6. An attempt 
to model the broadband X-ray spectrum from 3 to 200 keV 
showed that the data could be fit about equally well by a 
variety of models - thermal Comptonization with a coronal 
temperature of 40 keV and an optical depth of 3.6, a dif- 
ferent Comptonization model with a temperature of 87 keV 
and an optical depth of 1.6, or a cutoff power law with an 
index of 1.45 and a cutoff energy of 160 keV, plus an ad- 
ditional blackbody component, but several of the relatively 
good (i.e. x 2 l v < 2) fits suggested a covering fraction for 
Compton reflection consistent with zero (Dove et al., 1998). 
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What the past analyses lack and these analyses provide is 
full broadband coverage from 3 to 1000 keV. The combina- 
tion of Ginga and OSSE leaves a gap from roughly 30 to 50 
keV where the reflection component most affects the broad- 
band spectrum while the use of RXTE alone fails to reach 
the energies where the spectrum cuts off most steeply at ~ 
400 keV. Here we present the results of joint RXTE/OSSE 
fits to the spectrum of Cygnus X-l taken in December of 
1998. These data present the first gap-free energy spectra of 
Cygnus X-l from 3 to 1000 keV. 



2 OBSERVATIONS 

We examine data from three OSSE viewing periods where 
there was simultaneous RXTE data. The color changes mea- 
sured in the PCA and in the ASM over the one week periods 
are very small (less than a few percent; see Figure la for the 
PCA colors), suggesting that the spectrum remains stable 
over these time periods. Also, the ASM summed band inten- 
sities (1.3-12.0 keV) show a correlation timescale of about 
20 days (see Figure lb; the 1996 data have been excluded 
so the measure will not be biased by the state transition 
that occured then). The OSSE fluxes in the highest energy 
bands are weak and only detectors 3 and 4 from OSSE were 
turned on. Hence the full OSSE viewing periods, rather than 
just the strictly simultaneous data, are used as the average 
spectra for the analysis. The data are extracted using the 
standard reduction procedure from the CGRO science cen- 
ter. 

The RXTE data are analyzed using the standard RXTE 
GOF criteria for spectral extraction. Data from only the top 
layers are used, and all five PCUs are added together. The 
HEXTE data are also extracted using the standard reduc- 
tion criteria. The data used from each instrument is pre- 
sented in Table 1. 



3 ANALYSIS 

The data are then fit with the EQPAIR model using XSPEC 
11.0 (Arnaud 1996). The model allows for thermal and 
non-thermal comptonization of soft photons, pair produc- 
tion, Compton reflection, and bremsstrahlung emission. The 
physics of the model are described in detail in Coppi (1998). 
It has previously been used to fit spectra from Cyg X- 
1 (Gierlinski et al. 1997; Frontera et al. 2001), GX 339-4 
(Nowak, Wilms & Dove 2002) and GRS 1915+105 (Mac- 
carone, Coppi & Taam 1999; Maccarone 2001; Zdziarski et 
al. 2001). One percent systematic errors are added to all 
channels. Many of the parameters are frozen. The inclina- 
tion angle is assumed to be 55 degrees, in agreement with the 
most recent measurements (Sowers et al, 1998). The com- 
pactness, i bb (£ bb = L bb OT /m e c 3 Rcor, where L bb is the lu- 
minosity of the seed blackbody photon distribution, or is 
the Thomson cross-section, m e is the electron mass, and 
Rcor is the radius of the emission region) of the seed photon 
distribution is assumed to be unity. The elemental abun- 
dances are assumed to be solar. The reflector's temperature 
is assumed to be 100 eV. The overall normalizations of the 
HEXTE and OSSE components are allowed to float, but 
are extremely well constrained since even small deviations 



(i.e. ~ 1%) from the best fit normalizations will produce ex- 
tremely statisically significant discontinuities in the spectra. 
The PCA data are fit from 3 to 20 keV, the HEXTE from 17 
to 190 keV and the OSSE from 100 to 1000 keV. Additional 
model components of neutral hydrogen column absorption 
and a gaussian component to fit the iron emission line are 
also included. Since RXTE is not sensitive to the energies 
where absorption by the interstellar medium is most impor- 
tant, the neutral hydrogen column is fixed to 5 x 10 21 in 
accordance with ASCA and BeppoSax results (Gierlinski et 
al. 1999; Frontera et al. 2001 ). The iron line is restricted to 
fall between 6.1 and 7.0 keV in energy and to be less than 
1.2 keV in width. The parameters that are allowed to vary 
freely are the temperature of the seed blackbody photon 
distribution (kT bb ), the thermal electron compactness (£th), 
the reflection fraction, R and the ionization parameter £, as 
well as an overall normalization. Additionally, the HEXTE 
and OSSE data sets are allowed to have their normalizations 
float freely since the dead time correction for the HEXTE 
is uncertain (especially with the standard RXTE GOF dis- 
tribution response matrices which we have used) as are the 
cross-calibrations between RXTE and OSSE. Additionally, 
since the data from the two satellites are not strictly simulta- 
neous there may be some physical difference in the luminosi- 
ties apart from calibration issues. Typically, the HEXTE and 
the OSSE are both normalized to about 70% of the PCA. 
We also note that these observations are typical for Cygnus 
X-l as the power spectrum shows breaks at ~ 0.1 Hz and 
~ 5 Hz, which are approximately the median values for the 
break frequencies (Pottschmidt et al., 2001). Since the spec- 
tral indices tend to change much less than the variability 
timescales (see e.g. Belloni & Hasinger 1990), and the two 
are well correlated, this represents strong evidence that the 
observations are typical. 

3.1 Pure thermal Comptonization model fits 

The results of the fitting are presented in Table 2. A typical 
spectral fit is plotted in Figure 2. The residuals to this fit 
are shown in Figure 3. Most of the parameters vary little 
from observation to observation, in accord with the result 
that the colors in the PCA bands vary little over this time 
span. The parameters that change the most are the iron line 
equivalent width, the ionization parameter for the reflection, 
and the and the iron line's physical width. Since these pa- 
rameters are all most sensitive to the region between 6 and 8 
keV (the ionization level parameter depending primarily on 
the energy of the iron absorption edge), where there are not 
sufficient independent energy bins to fit the required number 
of parameters, one cannot be confident that the the differ- 
ences are real. Furthermore, the best fit equivalent widths 
of the iron lines are of order 100 eV for these observations. 
Previous measurements of the iron line for Cygnus X-l in 
the hard state from ASCA show ~ 20 cV equivalent width 
(Ebisawa et al. 1996). RXTE fits to the Crab spectrum can 
suggest the presence of an iron line with an equivalent width 
of about 80 eV, so features this weak should not be trusted 
regardless of the formal statistical significance of the result. 
Errors in the iron line measurements can lead to significant 
errors in the edge position. That the line energies are found 
to be below 6.4 keV at a strong signficance level (in dis- 
agreement with results from ASCA) also indicates that the 
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Figure 1. (a)The PCA colors over the first week of the observation. The colors change very little within the observations and by only a 
few percent from observation to observation. A few points are omitted from the plot near zero orbital phase where nH variations make 
a large difference in the colors and count rates (b) The autocorrelation function of the summed band (1.3-12.0 keV) ASM lightcurve. 



narrow feature is more likely related to errors in the response 
matrix than to a real physical feature. 

Other than the line parameters, a clear picture for the 
spectrum emerges where the disk temperature is fit to be be- 
tween 150 and 250 eV, the ratio of thermal compactness to 
seed photon compactness is about 9, and the optical depth is 
about 1.3. Thus previous results that the corona of Cygnus 
X-l in the hard state can be well fit as a photon-starved 
thermal plasma of moderate optical depth are verified. We 
have done several simulations and found that changes in the 
blackbody temperature can be measured due to changes in 
the curvature of the "power law" component of the spectrum 
at energies well in excess of the blackbody peak. For these 
observations, the blackbody temperatures vary very little; 
this model would be sensitive to variations about twice as 
large. Since slightly different models for the blackbody seed 
photon distribution can change the observed temperature 
and because of uncertainties in the RXTE response matrix 
and the lack of an RXTE response below 2 keV where the 
blackbody peaks, absolute measurements of the seed pho- 
ton temperature should be treated with skepticism. Relative 
changes, however, can be reliably measured. 

3.2 Hybrid Thermal/Non-thermal 
Comptonization model fits 

Additional fits are made with the inclusion of a non-thermal 
component to the spectrum. That is to say, the electron 
distribution then consists of a Maxwellian distribution (the 
"thermal" Comptonization component) and a power law 
component (the "non-thermal" Comptonization component, 
which is restricted to have electron Lorentz factors from 1.3 
to 1000 and has a compactness £ nt h and an electron spec- 



tral index of T such that - -E~ r ). The value of T is 
required to be greater than 1, the hardest electron spectrum 
that can be produced by standard acceleration processes. 
The best fit parameters are shown in Table 3. These fits 
are marginally better than those without the non-thermal 
component, with improvements in likelihood of factors of 
about 5 for VP 705 and VP 706 (which are best fit by a 
nearly purely non-thermal Comptonization model), and no 
improvement in VP 707. Since it is unlikely that a major 
qualitative change (i.e. whether the corona is dominated by 
thermal or non-thermal electrons) has occurred while the 
spectra look almost identical, we also attempt to fit a model 
to the data with a purely non-thermal electron for VP 707. 
We find that the data are consistent with this model as 
well (x /v — 1.04), with Fi n j—3.3, £ nt h~8.5, and the other 
paraemters consistent with their values for both the pure 
thermal and the hybrid model cases. 

The difficulty of distinguishing a thermal plasma from a 
hybrid one has been discussed previously. It has been noted 
that a hybrid plasma and a purely thermal plasma will pro- 
duce nearly identical spectra given the same optical depth 
and the same mean energy shift per scattering (Zdziarski, 
Coppi, & Lamb 1990; Ghisellini, Haardt, & Fabian 1993). 
Still theoretical models do predict differences between the 
two plasmas above ~ 400 keV (see Figure 4), since the purely 
thermal Comptonization models will show a nearly expo- 
nential decay, while the hybrid models will show a power 
law decay. Distinguishing between these two best fits would 
require sensitivity at the 4 x 10 -6 photons/cm 2 /sec level, 
which is beyond the capability of the existing OSSE obser- 
vations. A 10-<7 detection of this feature should be feasible 
with INTEGRAL in ~ 10 6 seconds, and with ASTRO-E2 
in about 10 5 seconds. Measurements from COMPTEL also 
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Figure 2. The unfolded thermal Comptonization fit to the data for OSSE viewing period 707, plotted in units. The data arc in 

open pentagons, the model fits in the solid line and the reflection component in the dashed line. The PCA component has been multiplied 
by 0.7 to make its relative normalization approximately the same as those for HEXTE and OSSE. 



seem to suggest the presence of an extended tail to the spec- 
trum (McConnell et al., 2000), but these observations arc 
not taken simultaneously with X-ray detections and also re- 
quire long integrations over which variability of the coronal 
temperature could be important. 



3.3 The Overall Compactness and the Size of the 
Emission Region 

We also investigate the range of overall compactnesses al- 
lowed by the data to determine whether constraints can be 
placed on the radius of the emission region and on the frac- 
tion of the coronal optical depth that is e~/e + pairs. We re- 
fit the model for VP 707 with a purely thermal electron dis- 
tribution, while varying £ bb . We find at the 90% confidence 
level that the total compactness of the system cannot ex- 
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Figure 3. The residuals to the fit for OSSE viewing period 707, plotted as A\ 2 versus energy. The PCA points are open pentagons, the 
HEXTE cluster points are open triangles, the HEXTE cluster 1 points are filled triangles, the OSSE cluster 3 points are open squares 
and the OSSE cluster 4 points arc filled squares. 



ceed ~ 300. This constraint requires that R CO r > 3 x 10 cm, 
which is not a strong constraint since this value is ~ IRsch 
for a 10 Mq black hole. We also find that the total compact- 
ness cannot drop below ~ .06 or the bremsstrahhmg cool- 
ing rate will dominate the Compton cooling rate, yielding 
a bremsstrahhmg emission spectrum which has too sharp 
a cutoff at high energies. This constraint requires that 
Rcor < 1.5 x 10 10 cm. This provides a spectral confirma- 
tion that the ~ few second hard time lags seen in variability 
analyses cannot be produced by light travel times in a large 



corona as had been suggested previously (see e.g. Kazanas, 
Hua, & Titarchuk 1997) and refuted on the basis of variabil- 
ity measurements (Maccarone, Coppi, & Poutanen 2000). 

We conduct the same analysis to find the minimum and 
maximum compactness for the purely non-thermal spectral 
fit to VP 707. The best fit is for £ bb = l-5±g;|. The tighter 
bound in the lower direction is because there is some intrin- 
sic curvature to the spectrum and some 7—7 pair production 
is needed to reproduce this. The tighter bound in the up- 
ward direction is because the annihilation peak at 511 keV 



© RAS, MNRAS 000, 1-8 



6 T. J. Maccarone & P. S. Coppi 




10 100 1000 

Energy (keV) 



Figure 4. The plot of the thermal and non-thermal fits to the spectrum from VP 706, folded through an idealized response matrix. 
Note that the differences are only apparent above ~ 400 keV. 

is much sharper for non-thermal Comptonization than for 
thermal Comptonization (i.e. in the thermal Comptoniza- 
tion case, the annihilation line is reddened and broadened 
significantly by Compton scattering off the thermal electrons 
- see Coppi 1998). These constraints require a corona with a 
size scale between 2 x 10 8 and 7 x 10 8 cm. This result is not 
very sensitive to letting y m in float freely, so if a purely non- 
thermal electron distribution can be proven to exist, fairly 
strong constrints can be placed on the size of the emission re- 
gion. A hybrid model (rather than purely non-thermal) pro- 



duces similarly tight constraints as the purely non-thermal 
model (4b = 2.2± 1 5 5 ). 

3.4 The reflection fraction 

The relatively low reflection fractions in Cyg X-l and other 
X-ray binaries have been explained as coming from a geo- 
metric effect where there is a "hole" inside the inner edge of 
the accretion disk from which no reflection comes (Gilfanov 
et al., 1999), a beaming effect due to a mildly relativistic 
outflow from the disk (Beloborodov 1999), smearing of the 
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reflected component as it passes back through the corona 
(Petrucci et al., 2001) or an underestimate of the reflection 
fraction due to a highly ionized disk (Young et al. 2001). 
We test the last of these mechanisms on VP 707 by allowing 
freezing the ionization parameter at a value of 10 5 . We find 
that even for this large ionization parameter, the reflection 
fraction remains ~ 0.35. Furthermore, the spectral fits are 
substantially poorer (A^ 2 ~ 150) than the best fits, which 
are obtained with low ionization parameters. We also ignore 
the PCA data from 5.0 to 9.5 keV and re-fit the model to en- 
sure that the iron line and edge have no impact on the best 
fit parameters and find that none of the parameters change. 
The reason for the discrepancy between this work and that 
of Young et al. (2000) is likely that the reflection fraction in 
the ASCA/Ginga data sets they fit depends largely on the 
location and strengths of the iron edges since ASCA has no 
sensitivity above 10 keV and Ginga has none above 30 keV, 
while in our fits, the reflection fraction depends largely on 
the curvature of the spectrum in the 30-50 keV range, since 
there are many more bins in this region than the 6-9 keV 
region where the iron line and edge are found. We have also 
fit the data after removing the OSSE data and have found 
essentially the same fit as for the PCA and HEXTE data 
alone. 



3.5 Robustness of the fits to changes in the frozen 
parameters 

We investigate the robustness of the spectral fits to changes 
of several parameters that have been left frozen. We find that 
changing the seed photon distribution to a diskPN model, 
which gives the sum of several blackbodies as from a pseudo- 
Newtonian disk (Gierlinski et al., 1999), rather than a pure 
blackbody yields a slightly higher optical depth (r ~ 1.4) 
and a higher thermal compactness (£ t h ~ 12) without mak- 
ing any qualitative changes to the spectrum or the overall 
goodness of fit. We find that for the nonthermal model, 7 m i n 
can be increased to ~ 3 without affecting the quality of the 
fit. For such high values of y m in, the electron injection dis- 
tribution approaches a 5-function, so these fits are probably 
unphysical, and it is unlikely that the non-thermal electron 
distribution is cut off much below 200 keV. 



4 CONCLUSIONS 

We find that the 3-1000 keV spectrum of Cygnus X-l in the 
hard state can be well fit by a purely thermal Comptoniza- 
tion model, but that a substantial (and perhaps dominant) 
non-thermal component cannot be ruled out. For a thermal 
model, we find typical temperatures of the corona to be ~ 
90 keV, with optical depths of ~ 1.3, and that for two of 
the observations, the data can be modelled slightly better 
with a purely nonthermal electron distribution. We find no 
evidence for the existence of an additional optically thick 
component in the spectrum. We find that the data are in- 
sufficient to differentiate between a pair-domianted corona 
and an electron/ion dominated corona. We find that the re- 
flection fractions are low (i.e. R/2tt < 1) because there is 
relatively little spectral curvature in the 30-50 keV range 
and not due to a ionization effects in the disk. Future obser- 
vations with instruments such as INTEGRAL and ASTRO- 



E2 should prove capable of breaking the degeneracy between 
the thermal and non-thermal electron distributions by mea- 
suring the emission at several MeV more accurately than 
current data sets which have provided simultaneous X-ray 
measurements. 
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Observat 
ObsID 

30158-01-01-00 
30158-01-02-00 
30158-01-03-00 
30158-01-04-00 
30158-01-05-00 
30158-01-06-00 
30158-01-07-00 
30158-01-08-00 
30158-01-09-00 
30158-01-10-00 
30158-01-11-00 
30158-01-12-00 
Viewing Period 

705 

706 

707 



ions Analyzed for 
Start Time 

10/12/97 07:07:53 
11/12/97 07:06:14 
14/12/97 08:48:14 
15/12/97 03:49:42 
15/12/97 05:26:14 
17/12/97 00:39:55 
20/12/97 07:11:18 
21/12/97 05:28:14 
24/12/97 23:03:14 
25/12/97 00:39:14 
30/12/97 02:18:03 
30/12/97 03:52:14 

Start Time 
1997-12-09 17:44:24 
1997-12-16 16:54:51 
1997-12-23 16:18:54 



This Work 

Stop Time 

10/12/97 08:30:14 
11/12/97 08:45:14 
14/12/97 10:20:14 
15/12/97 05:26:14 
15/12/97 07:09:14 
17/12/97 02:05:14 
20/12/97 08:29:14 
21/12/97 07:05:14 
24/12/97 00:39:14 
25/12/97 01:45:14 
30/12/97 03:52:14 
30/12/97 05:30:14 

Stop Time 
1997-12-16 14:08:23 
1997-12-23 13:27:23 
1997-12-30 14:15:46 



Table 1. The RXTE and OSSE observation times. VP 705 in OSSE coincides with RXTE observations 30158-01-01-00 through 30158-01- 
05-00, VP 706 coincides with 30158-01-06-00 through 30158-01-08-00 and VP 707 coincides with 30158-01-09-00 through 30158-01-12-00. 



Best Fit Parameters for the Thermal Model 

VP kT bb £ th t R ( F LE LW LN 

705 151±|§ 9.15±g| 1.26±-g| -23±-g^ 590±J™ 3.5xl0 37 6.1±£ -29±™ .0054±°°|} 

706 169±i| 9.12± : |g 1.31± : g§ .23±-°{ 486±ig° 3.7 xlO 37 6.1±;i -34± : §2 .0063± : ggii 

707 233±i| 8.49± : |i 1.30± : °§ -227± : |]^ 0.22± 3 £ 3.3 x 10 37 e.li, 1 , .97± : ^ .015±'° 3 



EqW 


X 2 /v 


60 


1.17 


70 


1.13 


200 


1.03 



Table 2. The best fit parameter values for the pure thermal Comptonization model. The error bars represent 90% confidence levels, 
while the double asterisks indicate parameters for which the 90% confidence level is outside the range allowed a priori. 



Zdziarski, A. A., Grove, J.E., Poutanen, J., Rao, A.R. & Vadawalc, 
S.V., 2001, ApJL, 554, 45 
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Best Fit Parameters for the Hybrid Model 



VP 


kT bb 


tth 


(■nth 


r 


T 


i? 




LE 


LW 


LN 


EqW 


x7" 


705 


192+J? 


0.18±-f| 




3-0±;i 


1.342±°°j 


9QQ + .003 

.ZOOZE 02Q 


41 0+ 190 


6-l±« 


■28±;i 


■ uuo ^ 3= .0014 


60 


1.12 


706 


197±f 


0.04± ; i| 


8.73±'J| 




1.367±°°5 


.ZOOIE 02 g 


16+J 00 


6-l±« 


i.o± ■} 


nl q+.003 


220 


1.07 


707 


258±§f 


6.38±i;i 




2.1±2, 5 


1 31+ - 02 


907+O.OI6 
•^ Ol=c 0.027 


5±f° 


6.1±2-2 


Q1+0.'23 


n i Q + .004 
.uidi 003 


180 


1.04 



Table 3. The best fit parameter values for the pure thermal Comptonization model. The error bars represent 90% confidence levels, 
while the double asterisks indicate parameters for which the 90% confidence level is outside the range allowed a priori. 
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